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The Coma cluster is one of the nearest galaxy clusters, and the first one in which 
a radio halo and a peripheral relic were discovered. While its halo and the central 
parts of the intracluster medium have been studied extensively, X-ray observations 
of the plasma near its relic have been scarce. Here, we present results from a re- 
analysis of a 22-ks archival XMM-Newton observation. Across the relic, we detect a 
shock of Mach number ^ 2. This excludes the previously suggested hypothesis that 
the relic was formed by turbulence. Furthermore, multiwavelenth observations and 
numerical models do not support the scenario in which the shock at the Coma relic is 
an outgoing cluster-merger shock. Instead, our results lend support to the idea that 
the relic coincides with an infall shock front formed just as the NGC 4839 group falls 
onto the cluster along a cosmic filament. 
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1 INTRODUCTION 

Over the last 40 years, diffuse, Mpc-scaled, steep-spectrum 
radio emission has been detected i n the outskirts of 
a number of galaxy c luster s (e.g., iFerrari et all l2008l : 
Ivan Weeren et al]|2009l . l201lh . These radio objects, named 
radio relics, have an irregular, arc-shaped morphology, gen- 
erally show a spectral index gradient across their surface, 
and are believed to be formed via acceleration at merger 
shock waves. 

The Coma cluster is one of the nearest galaxy clusters 
{z = 0.0231), and th e first one in which a peripheral radio 
relic was discovered (|BaIlarati et al.|[l98lh . However, X-ray 
observations of the plasma near the relic have been scarce. 
The Coma reli c, radio source 1253 -I- 2 75, has a length of 
roug hly 2 Mpc (jBrown fc RudnickllioTll ) . and is located SW 
of the cluster, at a projected radius of approximately 0.8r2ou. 
Using galaxy optical catalogues and 21-cm line observations 
take n with the 305-m r adio telescope of the Arecibo observa- 
tory, iFontaneUl (|l984h discovered that the galaxy distribu- 
tion in the region of the Coma cluster delineates a NE-SW- 
oriented cosmic filament, which includes the Coma cluster 
itself. In the context of matter accreting along this filament 
onto the Coma cluster, and given our current understanding 
of radio relic formation, the presence of a radio relic at the 
SW edge of Coma is not surprising. Yet, the formation mech- 
anism of the relic remains controversial |Ensslin et al.lll99j : 
iFeretti fc NeumannllioO^ : iBrown fc Rudnickll201lh . and un- 
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derstanding it requires multiwavelength, high-quality obser- 
vations. 

There has been only one attempt of detecting a shock 
at the Coma relic, using an XMM-Newton o bservation car- 
ried out in 2003 (jFeretti fc NeumannI [20061 ') . The analysis 
of the XMM-Newton data revealed a constant gas temper- 
ature of ~ 3 keV at and on the inner side of the relic. If 
the Coma relic can be explained by a shock front, then 
the temperature would be expected to increase at the relic, 
and then decrease rapidly by a factor of approximately 
(5M^ -\- 14)/16, where A4 is the shock Mach number (e.g 



iPaul et al] 20111). In the absence of a temperature increase 



or mmp. IFeretti fc Ne umann (200Q) suggested that the relic 
was instead formed by turbulence triggered as matter ac- 
cretes onto Coma. 

At the relic, the main source of X-ray emission is not the 
Coma cluster itself, but a smaller group of galaxies located 
along the NE-SW filament, NE of the relic . By comparing 
optica l and X-ray observations at the relic, [Neumann et al.l 
(|200ll ) have shown that the galaxies belonging to the group, 
including its brightest central galaxy, NGC 4839, are on av- 
erage situated approximately 300 kpc closer to the Coma 
centre than the group's X-ray halo. This displacement sug- 
gests that the s ubgro up is falling into the Coma cluster 

[Neumann et all (|200ll ). 

More recently, [Brown fc Rudnic3 ([201l[ ') studied the 
galaxy surface density and redshift distribution around the 
relic, using Sloan Digital Sky Survey (SDSS) data. They 
identified a 2-Mpc-long "wall" of galaxies on the inner side 
of the relic, infalling towards Coma's centre with a veloc- 



2 G. A. Ogrean et al. 



ity of approximately 500 kms^^. The length and position 
angle of the "wall" are very similar to those of the relic. 
They argue that the Coma relic traces a stationary infall 
shock formed behind the infalling group of galaxies. Infall 
shocks can form only in the early stages of a merger event, 
unlike (also stationary) accretion shocks, which occur due 
to continuous accretion onto the cluster. 

Here, we present results from a re-analysis of the XMM- 
Newton observation at the relic, which uses a more careful 
backg round treatment then the one of iFeretti fc NeumannI 
(|2006t) . Throughout the Letter, we assume a flat ACDM 
cosmology, with Ho = 70 km s~^ Mpc~^, Qm = 0.3, and 
f^A = 0.7. As X-ray emission at the relic comes mostly from 
the infalling subgroup, our analysis assumes a gas redshift 
of 0.0 246, equal to the redshift of NGC 4839 (|Smith et al.l 
which is the brightest cluster galaxy in the infalling 
group. At this distance, 1 arcmin translates to 29.7 kpc. All 
errors are la statistical errors, unless otherwise noted. 



Table 2. Best-fit background model. Temperatures are in units 
of keV. Normalizations are in default XSpec units per square 
arcmin. Power-law normalizations are measured at 1 keV. 

SPECTRAL MODEL: 

CALIB X SOLID_ANGLE X ( APEC + ( APEC + POW ) X PHABS ) 
Component F T Normalization 

LHB - 0.08+ (1.4 ± 0.23) X 10-'5 

GH - 0.18tro°oS S.9tHtxlO-^ 

CXB 1.4lt - 7.3j;o 54 X lO"'^ 

MOSl pn 

CALIB l-25±'(i;i; L02l°;°g 

f frozen 



2 OBSERVATIONS AND DATA REDUCTION 

This Letter is based on a 22-ks archival XMM-Newton ob- 
servation (ObsID 0058940701) centred on 1253+275 - the 
SW relic in the outskirts of the Coma cluster. The obser- 
vation was taken on June 10, 2003, using the thin filter. 
ICM emission from the Coma and the infalling subgroup is 
present across the whole field of view (FOV). Therefore, an 
additional observation was required to model the sky back- 
ground. We used ObsID 0406610301, a 15-ks, thin-filter ob- 
servation of the quasar HS 1251-1-2636, with an aim point 
angular distance of 0.92 degrees from the Coma relic, in 
the SW direction. Both datasets were reduced with the Ex- 
tended Source Analysis Software (esas) integrated in the 
Science Analysis System (SAS) version 12.0.1, and the most 
recent calibration files as of November 2012. The data pro- 
cessing steps closely follow those suggested in the XMM- 
Newton ESAS cookbook 0. 

Additionally, we used NASA's HEASARC X-ray Back- 
ground Tool to extract a ROSAT AU-Sky Survey (RASS) 
background spectrum from the same region as the spectrum 
of ObsID 0406610301. The RASS spectrum allows us to bet- 
ter constrain the LHB normalization. However, we note that 
excluding the RASS spectrum from the spectral fit does not 
siginificantly change any of the best-fit background parame- 
ters, with the exception of the LHB normalization; the LHB 
normalization contributes little to the fit though, given the 
very low LHB temperature. 

We checked for residual soft proton (SP) contamination 
by comparing the hard-band (10-12 keV for MOS, 12 - 14 
keV for pn) count rates inside and outside the field-of-view 
(FOV). iDe Luca fc Molendil (|2004h found that count rate 
ratios (in/out; TV) below 1.15 indicate event files that are 
essentially not contaminated by residual soft protons, while 
ratios between 1.15 and 1.3 translate to slight SP contami- 
nation. The higher the ratio, the larger the residual SP con- 
tamination. In Table [T] we list TZ for the fiare-filtered EPIC 
event files of the two observations, along with the "clean" 
exposure times. 

^ ftp : // xmm.esac .esa. int /pub /xmm-esas /xmm-esas . pdf 

^ http://heasarc.gsfc.nasa.gov/egi-bin/Tools/xraybg/xraybg.pl 



The quiescent particle background (QPB) was modelled 
spatially and spectrally, and subtracted from all the images 
and spectra. The modelling was done in three steps: (1) the 
spectral parameters of the detector corners are determined; 
(2) filter- wheel closed (FWC) observations with similar spec- 
tral parameters are identified in the archived-observation 
database; (3) the FWC data is scaled by the ratio of the 
observation corner spectra to FWC corner spectra. Finally, 
ESAS creates images and spectra of the QPB using the scaled 
FWC data. 



3 ANALYSIS 

3.1 The sky background 

Sky background spectra were extracted from the ObsID 
0406610301 event files in circular annuli between 6 and 12 
arcmin around the optical axis, and binned to a minimum of 
30 counts per spectral channel. The spectra were modelled 
as the sum of Local Hot Bubble (LHB), Galactic Halo (GH), 
Cosmic X-ray Background (CXB), and residual SP (RESP) 
emission. The LHB and the GH were described by thermal 
components, while the CXB and the RESP were described 
as power-laws. For the absorbed components (GH, CXB), 
the X-ray column density was fixed to the weighted average 
Galactic Hi column density in the Leiden/ Argentine/Bonn 
(LAB) Survey. Two constants were introduced to describe 
the solid angle corresponding to the extraction region, and 
to account for calibration errors between the EPIC detec- 
tors (calculated with respect to M0S2). Th e spectral model 
is similar to that of ISnowden et al] (|2008l l. and described 
there in more detail. Table [2] summarizes the best-fit sky 
background model. Despite our attempt to model RESP 
emission, none was detected in any of the spectra; this is 
not necessarily surprising, given that, within the error bars, 
TZ of all detectors is consistent with < 1.15. The fit was ob- 
tained with XSpec v12.7.1 in the energy band 0.5 — 7 keV 
(excluding the 1.2 — 1.9 keV band, which is contaminated 
by fiuorescent instrumental lines), and has x^/d.o.f. = 1.02. 



X-ray evidence of a shock at the Coma relic 3 



Table 1. Level of residual SP contamination and "clean" exposure times. TZ describes the ratio of count rates inside and outside the 
FOV, in the hard-energy band. See text for details. 



ObsID 


n MOSl 


7^ MOS2 


It pn 


Exp. time MOSl 


Exp. time MOS2 


Exp. time pn 


0058940701 


1.272 ±0.094 


1.046 ± 0.070 


1.207 ±0.119 


16.0 


16.5 


9.2 


0406610301 


1.074 ±0.081 


1.167 ±0.081 


1.210 ±0.115 


10.3 


10.3 


6.3 



3.2 Evidence of a shock 

To the best-fit background model above, we added an ab- 
sorbed, single-temperature thermal component to model 
emission from the ICM. Additionally, we used a power-law 
to model RESP emission. The X-ray column density was 
fixed to 8.06 x 10^® atoms cm~^ - the weighted average 
LAB value measured in a circle of radius 1 degree around 
the relic. The metallicity was frozen to 0.2 solar, while the 
redshift was set to the redshift of NGC 4839, z = 0.0246 
l|Smith et al.ll200GD . The factor that accounts for calibration 
errors was calculated with respect to the MOS2 spectrum. 
The RESP power indices and normalizations were coupled 
between spectra of the same detector, taking into account 
that the normalizations differ only by a constant factor that 
can be calculated from the output of the xmm-esas routines. 

The count statistics of the XMM-Newton observation 
only allow measuring the ICM temperature to within ^ 10 
percent in about four regions. We chose them to cover a re- 
gion on the relic, and regions in front and behind the relic. 
Because the surface brightness appears enhanced to the SE, 
we modelled the ICM in this region separately. The regions 
are shown in Figure (2] The spectra were all grouped to a 
minimum of 30 counts per spectral channel, and fitted simul- 
taneously in the same energy band as the sky background 
spectra. By fitting all the spectra simultaneously, the RESP 
parameters, the only ones linked between spectra, are better 
constrained. The best-fit spectral model for each region is 
summarized in Table[3] No SP component has been detected 
in the pn spectra [IZ for the pn data is consistent with 1.15 
within the error bars). The best-fit RESP parameters for the 
MOSl spectrum of the inner 12 arcmin around the optical 
axis are F = LltHf and M = 3.7t^j^ x 10"^ with the 
normalization in default XSpec units. The best-fit calibra- 
tion factors ranged between 0.73 and 0.96. The overall fit 
has xVd-o.f. = 1.00. 

Figure [T] shows the fitted spectra of the four regions. 
The NE and central regions have roughly equal temperature, 
ss 3 keV , consistent with the previ ous temperatures calcu- 
lated by iFeretti fc Neumamil (|2006t ) for regions NE of and 
overlapping the relic. However, we measure a strong tem- 
perature jump at the SW edge of the relic; the temperature 
here decreases to « 1 keV, which points towards a shock of 
Mach number 1.81q'32, derived using the Rankine-Hugoniot 
relations for 7 = 5/3. 

We verified the robustness of our results by varying each 
background parameter within la of its best-fit value, and 
refitting the ICM spectra with the new background model. 
The newly fitted temperatures are all consistent with the 
previous best-fit values, within their Icr errors. Further- 
more, using data from the Solar Wind Ion Composition 
Spectrometer on board of the Advanced Composition Ex- 
plorer satellite, we also plotted the O^^/O"''^ ion ratio as a 



Table 3. Best-fit ICM models. Temperatures have units of keV, 
while normalizations (A/") are given in units of cm~^ arcmin"'^. 
For compari son, we also list the best- fit temperatures (T™^^) 
obtained by iFeretti fc NeumannI 1 I2OO6I ) for regions on the inner 
side of the relic (NE), and on the relic (centre); however, while 
their regions are included in ours, the exact region selection is 
different. 

SPECTRAL MODEL: 

CALIB X SOLm_ANGLE X ( PHABS X APEC ) 



Region T AT T^Noe 



NE 


2.9 ±0.21 


(1.7 ±0.073) X 10 


-5 


o O+0.7 
■^"^-0.5 


centre 


2.9 ± 0.21 


(9.3 ±0.36) X 10- 


6 


0+0.3 

0-0„g 4 


SW 


1 fi+0.20 
^^■"-0.32 


(2.0 ±0.30) X 10" 


6 




SE 


4+0.42 


(1.3 ±0.075) X 10 


~5 




12 arcmin 


9 fi+0.12 
^■"-0.11 


(8.2 ±0.21) X 10" 


6 






20 30 40 

r (arcmin) 



Figure 3. Background-subtracted, vignetting-corrected surface 
brightness profile across the relic. The surface brightness was cal- 
culated in the energy band 0.4—4 keV. Error bars show Poissonian 
errors. Dotted lines show the projected extent of the relic. The 
dashed horizontal line represents the background level (sky plus 
instrumental). Overplotted are best-fit models corresponding to 
an underlying power-law density profile (red; x'^/d.o.f. = 7.9), 
and two power-laws with equal indices and a jump fixed at a 
radius of 35 arcmin (blue; x^/d.o.f. = 3.0). 



function of time for one week before and after the date of 
the Coma and background observations, to see whether the 
data was affected by solar wind charge exchange (SWCX) 
(jSnowden et al.ll2004h . The 0+''/0+® ratio during both of 
the XMM-Newton observations was among the lowest in the 
two-week period, with values < 0.1. 
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NE REGION (1.09) CENTRAL REGION (0.95) SW REGION (0.90) SE REGION (1.06) 




Energy (keV) Energy (keV) Energy (keV) Energy (keV) 



Figure 1. Top: MOSl (black), MOS2 (red), and pn (green) spectra of the NE, central, SW, and SE region, with overplotted best-fit 
models. Bottom panel shows the ratio of the data to the model. The power-law corresponds to the MOSl RESP component, which was 
added to the model without being folded through the instrumental effective area. Reduced values for each fit are given in the plot 
title, in parenthesis. 
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Figure 2. Left: Instrumental background-subtracted surface brightness map of the ICM near the Coma relic, binned by a factor of 2 and 
adaptivcly smoothed. The image corresponds to the energy band 0.4 — 4 keV. Overlaid are WSRT 90-cm radio contours. The circular 
region has a radius of 12 arcmin. Centre: The same X-ray image, now showing the partial annuli used for the surface brightness profile. 
Right: Temperature map. Black lines separate the regions used for spectral analysis. 



3.3 No surface brightness jump 

AH shocks associated with radio relics were first iden- 
tified based on discontinuitie s in t he surf ace brightness 
(e.g., A3667. lFinoguenov et al.1 (|2010h : A754. iMacario et al.1 
()201l| )l. We extracted a surface brightness profile across 
the relic in partial annuli with an angular opening of 15 
degrees and with a width of 15 arcmin, centred on NGC 
4839. The profile is shown in Figure O It does not show any 
clear density discontinuity at the relic. Under the assump- 
tion of spherical symmetry, however, the surface brigthness 
profile is not consistent with a simple power-law density pro- 
file within the FOV (because the cluster's redshift is low, 
there is essentially no difference between a beta profile and 
a power-law profile on length scales equal to the XMM- 
Newton's 30 arcmin FOV). Describing instead the density 
profile by two power-laws with a jump at a certain shock ra- 
dius yielded improved, but still unsatisfactory fits. We tried 
fixing the density jump to the shock compression expected 
from the temperature-derived Mach number, but the fit re- 
mained poor. The best fit was obtained for power-laws with 
equal indices (F — 1.07), and a jump fixed at 35 arcmin; the 
shock compression ratio was approximately 1.2, but the fit 
had a x^/d.o.f. of only 3. 

In conclusion, the surface brightness profile does not 



appear consistent with ~ 2 or with an underlying beta- 
model density profile, at least not for the simple spherical 
models assumed here. 



4 THE NATURE OF THE SHOCK 

5 W of the Coma relic, we have identified a temperature jump 
that corresponds to a Mach number of LSlg ^j. However, no 
similar discontinuity is found in density, assuming a simple 
spherical-shock model. This result is not necessarily surpris- 
ing though. At the position of the relic, the observed surface 
brightness profile is the integrated emissivity of both Coma's 
ICM, and the ICM of the infalling NGC 4839 group. This is 
a first deviation from our simplistic model, which will flatten 
the observed density profile. Furthermore, if the infall does 
not occur in the plane of the sky, then projection effects will 
smooth out the apparent density jump. Further smoothing 
can be caused by substructure along the line of sight. A 
similar effect has been seen in the merger A665, in which a 
shock region in front of the moving subcluster core has been 
identified in the temperature map, but is not associated with 
a surface brightness disc ontinuity (|Markevitch fc VikhlinirJ 
I2OOII : iGovoni et aLll2004l ) . 

The temperature jump measured at the relic corre- 
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spends to a shock of A4 2. The nature of the shock/rehc 
has been longly debated. lEnssHn et alj (|l998h suggested that 
the rehc traces an accret ion shock caused by ongoin g gas ac- 
cretion onto the cluster. iFeretti fc NeumannI (|2006l ) did not 
detect any evidence of a shock at the rehc, and suggested 
that it might hav e been formed by turbulence. Recently, 
I Brown fc Rudnic^ (|201lh proposed an alternative explana- 
tion: that the Coma relic traces an infall shock, the only 
infall shock claimed so far, caused by the infall of the NGC 
4839 group onto Coma. Another possibility is that the shock 
is another more common outgoing merger shock. 

It is now clear from optical and X-ray data that the 
NGC 4839 group is falling into Coma along a NE-SW- 
oriented cosmic filament, at a velocity of ~ 1500 kms^^ 
(|Colless fc Dunnlll996l : iNeumann et al.ll200ll ). The distance 
between the relic and NGC 4839 is approximately 900 kpc, 
which means that the group crossed the current position of 
the relic approximately 0.5 Gyr ago. If the shock is an out- 
going merger shock, the temperatures on both sides of the 
relic imply a shock velocity of roughly 1200 kms~^. In this 
geometry, 0.5 Gyr ago the infalling group was at the current 
position of the relic, while the shock was in-between Coma 
and the NGC 4839 group, 0.6 Mpc ahead of the group. This 
cluster/group configuration would be consistent with the ini- 
tial stage of a binary galaxy cluster merger predicted by nu- 
merical simulations, when a shock is formed at the interface 
between the two clusters. Numerical simulations also show 
that, following first core passage, outgoing merger shocks 
will be released into the ICM, outward of the two cluster 
cores; outgoing merger shocks are not produced behind the 
infalling subcluster core during its first infall onto the more 
massive cluster (e.g..lTakizawall2000l : lRicker fc Sarazi nll200ll : 
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The NGC 4839 group is currently 
located approximately 1.5 Mpc away from Coma's centre, 
and it is observed well before first core passage. Therefore, 
the shock at the Coma relic, observed behind the infalling 
group, cannot be an outgoing merger shock triggered during 
a simple binary merger event between the group and Coma. 

The low Mach number of the shock also excludes the 
possibility that we are observing an accretion shock pro- 
jected at a radius smaller than Coma's virial radius. Accre- 
tion shocks are exp ected to have Mach numbers > 10 (e.g., 
iMiniati et al.ll200ol '). while our XMM-Newton data only re- 
veals a shock of ~ 2. 

The Coma relic ha s an integrated spectral index of 1.18 
(|Thierbach et al.ll2003l '). Taking into account electron aging, 
the Mach number predicted by the injection spectral index 
is 3.5. The temperature discontinuity points towards a shock 
of Mach number ~ 2. However, colder gas from the cosmic 
filament must have been dragged into Coma's ICM in the 
wake of the infalling group, while the ICM of the group 
must have been ram-pressured stripped. Therefore, our as- 
sumption that the pre-shock and post-shock regions are best 
described by single-temperature thermal components is sim- 
plistic, and the actual strength of the shock is likely higher. 
Unfortunately, the low count statistics of our data do not al- 
low testing two-temperature spectral models. Detailed spec- 
tral index and polarization mapping of the relic, as well as 
numerical simulations aimed at modelling the group's infall, 
are needed to further characterize the shock at the Coma 
relic. However, by finding evidence that the relic traces a 
shock front, we can confidently eliminate the scenario in 



which the relic was formed by turbulence. Of the remain- 
ing possibilities, the most likely one appears to be that the 
shock at the Coma relic is the only infall shock detected so 
far. 
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